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ABSTRACT 

The results of a deep submillimetre observation using SCUBA of the powerful radio 
galaxy 3C324, at redshift z — 1.206, are presented. At 850 microns, emission from 
the location of the host radio galaxy is marginally detected at the 4.2a level, 3.01 ± 
0.72 mJy, but there is no detection of emission at 450 microns to a 3cr limit of 21 mJy. 
A new 32 GHz radio observation using the Effelsberg 100 m telescope confirms that the 
sub-millimetre signal is not associated with synchrotron emission. These observations 
indicate that both the mass of warm dust within 3C324, and the star formation rate, 
lie up to an order of magnitude below the values recently determined for radio galaxies 
at z ~ 3 — 4. The results are compared with dust masses and star formation rates 
derived in other ways for 3C324. 
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1 INTRODUCTION 

Searches for emission from distant galaxies at millimetre and 
submillimetre wavelengths provide a direct method of inves- 
tigating the star formation activity of these galaxies, since 
a large proportion of the optical and ultraviolet luminos- 
ity is reprocessed by dust to these longer wavelengths (e.g. 



dio sources with redshifts z ~ 1 are well-formed ellipticals 
whose stellar mass is dominated by an old stellar population 
(e.g. Stockton et al. 1995), and which have settled down into 



a stable r 1 ' 4 law profile (Best et al. 1998b). These galaxies 
are amongst the most massive at their epoch, with stellar 
masses of a few times 1O 11 M0, and oft en appear to lie a t the 



Soifer & Neugebauer lyyl). In addition, the strong negative et~al 



centre of you ng or forming clusters (Dickinson 1997 



K-correcliou caused by Lhe steepness uf the Ray leigh- Jeans 

tail of the thermal dust spectrum means that at submillime- 
tre wavelengths equivalent objects will appear as bright a t 
redshift z ~ 10 as they do at z ~ 1 (Blain & Longair 1993). 



Best 



1998b). Although a proportion of the aligned emis- 



Powerful radio galaxies can be observed out to extreme red- 
shifts, z > 4 (e.g. Rawlings et al. 1996), and therefore pro- 
vide a probe through which the star formation history of 
massive galaxies can be investigated. 

The optical and ultraviolet morphologies of powerful 
radio galaxies with redshifts z > 0.6 are dominated by emis- 
sion elongated and aligned along the direction of the radio 
axis ( |McCarthy et al. 1987| ; phambers et al. 1987[ ). The pre- 
cise details of this so-called 'alignment effect' remain un- 
clear, with scattered quasar light, nebular continuum emis- 
sion and jet-induced star formation all likely to contribute 
at some level (e.g. see Rottgering & Miley 1996). Underly- 
ing this aligned emission, however, the host galaxies of ra- 



sion may be associated with young stars, forming at a rate 
of order 100M Q yr" 1 (e.g. Best et al. 1996, 1997a, Cimatti 
et al. 1996) , the few galaxies that have been studied in detail 
so far show no strong signatures in the rest-frame ultravi- 
olet / optical spectra that would indicate star-formation at 
significantly higher rates (Cimatti et al. 1996, 1997). 

At redshifts z > 2, powerful radio galaxies have much 
more clumpy optical morphologies, with extreme examples 
showing many separate brig ht emission regions wit hin a spa- 
tial extent of 50 to 100 kpc (IPentericci et al. 19981). The in- 



dividual emission clumps have sizes of 2 to 10 kpc, and their 
profiles and colours are similar to those of the UV-dropout 
galaxies at z ~ 3 ( Giavalisco et al. 1996 ). Although these 
radio galaxies overall display a strong alignment effect, a sig- 
nificant proportion of the individual clumps lie away from 
the radio axis, suggesting that their blue colours are not di- 
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rectly associated with the alignment effect, but that these 
are more likely to be star-forming objects, each with star 
formation rates of order 10 solar masses per year. 

At redshifts z > 3, targeted millimetre and submillime- 
tre observations have succeeded in detecting a number of 
powerful radio galaxies (see Hughes et al. 1997 for a review, 
Ivison et al. 1998, Cimatti et al. 1998, Rottgering et al. 
1998, Hughes & Dunlop 1998). The masses of warm dust 
implied by these observations are a few times 10 8 Mq and 
the corresponding star formation rates are extreme. If it is 
assumed that all of the dust is heated by young stars, then 
thousands of solar masses per year of star formation must 
be on-going. These results are frequently interpreted as be- 
ing the burst of star formation which forms the bulk of the 
stellar mass of these galaxies. Supporting evidence for this 
hypothesis comes from a deep Keck spectrum of the radio 
galaxy 4C41.17 (z = 3.8) at rest-frame ultraviolet wave- 
lengths, which shows strong absorption features indicating 
that the ultraviolet continuum of this galaxy is dominated 
by hot young sta rs forming at rates of up to llOOM© yr _1 
( |Dey et al. 1997| ). 

Despite the successes at the highest redshifts, submil- 
limetre observations of powerful radio galaxies at redshifts 
2 ~ 1 have been sparse and unsuccessful (e.g. see Hughes 
& Dunlop 1998). Studies of these galaxies are of great im- 
portance, firstly for understanding the continued evolution 
of the interstellar medium of massive galaxies, and sec- 
ondly for investigating the extent to which the powerful ac- 
tive galactic nucleus (AGN) can be responsible for heating 
the dust within the radio galaxies. The advent of highly 
sensitive instruments such as the Submillimetre Common- 
User Bolometer Array (SCUBA; Holland et al. 1998) on the 
James Clerk Maxwell Telescope^] (JCMT) has opened up 
the possibility of more detailed studies of these objects. 

In this letter we present submillimetre observations of 
one such radio galaxy, 3C324 at redshift z = 1.206. The 
very strong alignment effect and the scattering properties of 
this galaxy, together with a central absorption feature in the 
HST image, make it one of the most promising candidates 
in t he 3CR z ~ 1 sample for a large dust mass (Best et al 
1996; Cimatti et al. 1996). In Section we describe the ob- 



servations of 3C324 and the data reduction. In Sectional 
present our results and discuss their implications. Through- 
out this letter we adopt qo=0.5 and Ho = 50kms~ Mpc~ . 



2 OBSERVATIONS 

2.1 SCUBA observations 

3C324 was observed using SCUBA on the JCMT for a to- 
tal of 105 minutes on source during the nights of November 
30th and December 2nd, 1997. SCUBA has two arrays of 
bolometric detectors, cooled to 0.1K. The short-wavelength 
(SW) array has 91 feedhorns and can be operated at 350 
or 450 /im; the diffraction-limited beamwidth at 450 /im has 
a half-power diameter of 7.5 arcsec. The long-wavelength 
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(LW) array has 37 feedhorns and can be operated at 750 or 
850 /im; at 850 /im the beamwidth is 14 arcsec. A dichroic 
beam-splitter enables observations using both arrays simul- 
taneously. 

In the standard SCUBA photometry mode, observa- 
tions are made using the central pixels of each array, which 
are aligned with each other to within about an arcsecond. A 
jiggling procedure is employed whereby the secondary mir- 
ror is jiggled around a 9-point filled square, with a 2 arcsec 
offset between each pointing. This method has been shown 
to produce the best photometric accuracy (cf Ivison et al. 
1998) . The integration time for each pointing of the jiggle is 1 
second. During this period the secondary mirror is chopped 
at 7 Hz between the source and a reference sky position, 
usually with a chop-throw of between 60 and 180 arcsec in 
azimuth. Following the 9 second jiggle, the telescope is nod- 
ded so that the chop position is placed at the opposite side 
of the source and the jiggle pattern is repeated. 

The observations of 3C324 were made at 850 /im and 
450 /mi using the two-bolometer photometry mode, which 
is an adaptation of the standard photometry mode. In this 
mode the chop-throw is not fixed in azimuth, but is fixed 
such that, in the primary nodding position, the source is 
chopped to the centre of one of the other bolometers of 
the LW array. This 'chop bolometer' is chosen to be one 
for which the chop-throw lies as nearly as possible in az- 
imuth, and whose position corresponds closely to that of a 
bolometer in the SW array. In two-bolometer mode, the 
source is thus observed for half of the total integration 
time in the central bolometer, plus a further quarter in the 
chop bolometer, improving the observing efficiency. Due to 
the apparent curvature of the arrays on the sky, the chop- 
throw in the reference nodding position does not centre the 
source on a bolometer. A three-bolometer photometry mode 
would maximise the observing efficiency and is planned for 
SCUBA, but was not operational at the time of these obser- 
vations. 

350 18-second integrations were made of 3C324, split 
into 7 sets of 50 integrations. Telescope pointing checks were 
made using 1611+343 before each set of integrations; the off- 
sets were small, typically below 2 arcsecs. Skydips were made 
at intervals of about 2 hours throughout both nights, and 
showed the sky to be stable. On November 30th the atmo- 
spheric zenith opacity was consistently about 0.12 at 850 fim 
all night and between 0.55 and 0.65 at 450 /mi. On Decem- 
ber 2nd the opacity was about 0.25 at 850 /mi and 1.6±0.3 
at 450 /tm. The observations were made whilst 3C324 was 
at low airmass (mean value 1.24). Calibration observations, 
using 3C273 on the first night and IRC+10216 on the sec- 
ond, were taken using the same two-bolometer photome- 
try mode, the fluxes of these objects being bootstrapped 
to Mars. The calibration uncertainty for these data is esti- 
mated to be < 10% at 850/im, but possibly as much as 25% 
at 450/im. 

Th e data were reduced usin g the SCUBA software, 
SURF (Jenness & Lightfoot 1998). The reference measure- 



ments were subtracted from the signal beams, the bolome- 
ters were flatfielded using the standard SCUBA flatfield, and 
the extinction correction was applied. Noisy integrations and 
strong spikes in individual bolometers were rejected; this re- 
moved a little under 5% of the data. The residual sky back- 
ground was removed by subtracting the mean signal from 
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the off-source bolometers in the inner ring around the cen- 
tral bolometer, excluding those bolometers which had sig- 
nificantly higher than average noise (see also Ivison et al. 
1998). 

The integrations were then concatenated to form a sin- 
gle dataset, and the consistency of this dataset was investi- 
gated using a two sample Kolmogorov-Smirnov (KS) test. 
In this way, periods when the changes in, for example, the 
atmospheric conditions or the telescope focus may hav e af- 
fected the data quality can be identified and removed (Jen- 



& iLightfoot 199^)7 The data was split into subsamples, 
and the first two samples were compared for consistency. 
The second sample was rejected if the probability that the 
two samples are drawn from the same parent sample was 
below a given limit; otherwise, it was concatenated with the 
first sample. The resulting sample was then compared with 
the third subsample, the process was repeated, and so on. 
For each wavelength the KS test was run on both the central 
bolometer and the chop bolometer with a variety of inputs: 
the number of subsamples, the order in which they were sup- 
plied, and the limit for rejecting the samples were all varied. 
The results obtained were consistent, and rejected about 5% 
of the data. 

At 850 fim, the signal measured in the central bolometer 
was 3.65 ± 1.17mJy, and that in the chop bolometer 4.46 ± 
2.55 mjy (the error in the latter measurement is higher since 
this bolometer was only on-source for half as long as the 
central bolometer). Combined, these give 3.78 ± 1.05 mjy, a 
signal with a 3.6<r significance. 3C324 has a lso been observed 
using SCUBA by Hughes & Dunlop ( [l998| ). They determine 
a flux density at 850 /im of 2.4 ± 1.0 mjy. Combining this 
with our data gives a combined signal of 3.01 ± 0.72 mjy, 
a marginal detection at the 4.2<r level. At 450 /J,m no signal 
was detected to a 3<r upper limit of 21 mjy. 

2.2 Effelsberg Observations 

3C324 was observed using the 100 m Effelsberg Telescope on 
July 15th 1998, at a central frequency of 32 GHz. The 3-feed 
receiver system installed in the secondary focus was used in 
a multi-beam mode. Each horn feeds a 2-channel receiver 
with an IF polarimeter providing full Stokes information 
simultaneously. The bandwidth was 2 GHz. 

The observations were made using a cross-scan in the 
equatorial coordinate system, with the main beam scanning 
a distance of 4 arcmin at a scanning speed of 10 arcsec per 
second. The offset feeds were used to efficiently remove at- 
mospheric noise. 32 such sub-scans were made of 3C324 (16 
in a north-south direction, and 16 east-west), and the data 
in the combined scan was sampled at 3.8 arcsec intervals. 

4 sub-scans of the source 3C286 were also made, and 
used to calibrate the flux density of 3C324, taking the flux 
density of 3C286 at 32 G Hz to be 1.779 Jy according to the 
scale of Ot t et al. (1994), which corresponds to the Baars 
scale (1977) at lower frequencies. The flux density of 3C324 
at 32 GHz was thus determined to be 32.4 ± 9mJy. 



3 DISCUSSION 

In Figure |l] we plot the results of these SCUBA observa- 
tions, together with the flux densities of the synchrotron 



emission at radio wavelengths from the literature (using the 
compilation of Herbig & Readhead 1992, and new 5 and 8 
GHz measurements from Best et al. 1998a), and our new 
Effelsberg data. 

The radio synchrotron emission of double radio sources 
steepens at high frequencies due to electron cooling, the 
precise details of the steepening depending upon the syn- 
chrotron ageing model adopted. The radio spectrum of 
3C324 was compared with the three popular age ing models, 
and the best fit was pro vided using a KP model ( Kardashev 



1962 



Pacholczyk 1970 ). This is shown on Figure [jj the fit- 



ted synchrotron spectrum passes more a factor of 10 below 
the 850 /im SCUBA signal. It is beyond the scope of this 
letter to discuss in detail the synchrotron ageing models: we 
merely note that although the KP model reflects a some- 
what unphysical situation, since it does not all ow a uniform 
distribution of pitch angles, Carilli et al. (1991) showed that 
this model also provides the best match to the lobe emission 
of Cygnus A, a low redshift analogue of 3C324; therefore, it 
provides an acceptable template. 

The hotspots and radio core, where the electrons are 
continuously injected, may, however, show less spectral 
steepening. The radio core flux densities at 5 GHz and 8 GHz 
are only 0.17 and 0.14 mjy respectively (Best et al. 1998a; 
see Figure |lj), and so its flux density at 850 /im will be neg- 
ligible. The hotspots, which lie towards the extremes of the 
SCUBA 850 micron beam (the radio source angular extent 
is 11.5 arcsec, compared to a beam FWHM of 14 arcsec), 
cannot be resolved from the lobe emission even in our high 
resolution 8 GHz maps, and so precisely determining their 
contribution to the radio source flux density is not possible. 
The poor agreement of a continuous injection synchrotron 
model with the 15 and 32 GHz data points, however, sug- 
gests that they do not dominate the flux density even at 
the highest radio frequencies. In conclusion, the shape of 
the high frequency radio spectrum appears to rule out the 
high frequency tail of the synchrotron emission as a possible 
origin of the 850 fim signal. 

Although we cannot categorically state that a signal 
with 4.2a significance constitutes a detection, the agreement 
between the signal we derived and that deriv ed fr om an in- 
dependently data set by Hughes and Dunlop ( 199J ) , and the 
consistency between the two bolometers and stability of the 
signal throughout the seven sets of 50 integrations in our 
observations (a positive signal was measured in ten of the 
fourteen measurements), adds support to its reality. For the 
remainder of this letter we shall perform calculations based 
upon an 850 /im flux density of 3.01 mjy: it should be re- 
alised, however, that the values derived should strictly be 
treated as upper limits instead of derived quantities. 

The mass of warm dust, Md, can be calculated from the 
sub-millimetre flux using the equation: 



S (fobs 



(1 + z)K d (v [cs t)B(v leat ,Td) 



where S is the observed flux density, u h s and f rcst are the 
observed and rest-frame frequencies, Dl is the luminosity 
distance, z is the redshift, red is the mass absorption coef- 
ficient, B is the black-body Planck function, and Td is the 
dust grain temperature. To allow c omparison with the dis- 
cussion of the z > 3 radio galaxies (Hughes et al. 1997), we 
adopt a dust temperature Td = 50 K, and a mass absorption 
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Figure 1. The spectral energy distribution of 3C324. The open diamonds represent the current SCUBA data, and the solid, dotted 
and dashed lines represent isothermal grey-body emission with an emissivity index (3 = 2 for dust at temperatures of 70K, 50K and 
30K respectively. The filled circles are the radio flux densities of 3C324, taken from Herbig & Readhead (1992) and Best et al. (1998a). 
The open triangle is our new 32 GHz data point. The dash-dot line shows a single power-law fit to the synchrotron emission at radio 
frequencies below 10 GHz. The open circles are the flux densities of the radio core emission only, taken from Best et al. (1998a). 



coefficient « d = 0.067(^ rost /250GHz)' 3 m 2 kg" 1 with (3 = 2. 
This provides a mass of warm dust in 3C324 of 1.2 x lO 8 A/0. 
If it is assumed that this dust is heated primarily by young 
stars, then a simple scaling between the star-formation rate 
and the submillimetre lumino sity can be obtain ed using 
nearby starbursts such as M82 ( [Hughes et al. 1994 ), provid- 
ing a current star formation rate for 3C324 of 35OM0yr _1 . 

Adoption of a lower temperature (the 450 /im SCUBA 
upper limit implies that if the 850 /im is real then the dust is 
at a temperature T < 45 K) would increase the dust mass, by 
about a factor of 2 for Tk = 30 K. The dust mass would also 
increase if a flatter grain emissivity index were used (about 
a 50% increase for (3 = 1.5). A full discussion of the uncer- 
tainties of these values can be found in Hughes et al. (1997). 
The du st masses derived are also strongly dependent upon 



Ho and go through the luminosity distance dependence: for 
Ho = lOOkrns - Mpc - the mass would be a factor of 4 
lower than that derived here; for go = it would be nearly a 
factor of 2 higher. All these conversion factors would, how- 
ever, apply similarly to any star formation rates determined 
from the optical and ultraviolet emission, and (in the case 
of qo, to an even greater extent) to the dust masses derived 
for the highest redshift radio galaxies. 

The derived mass of warm dust within 3C324, 1.2 x 
10 s Mq, can be compared to dust masses derived in other 
ways for this radio galaxy. HST images of the galaxy show 
bright extended ultraviolet emission, but the central re- 
gions of the galaxy are ob scured, probably by a dust lane 



with E(B - V) > 0.3 (Longair et al. 1995 



Dickinson 



can be related to the extinction using the equation M<j = 
T,{Ab) /Fb, where E is the area covered by the dust extinc- 
tion, (Ab) « 4E(B — V) is the mean B-band extinction, and 
Tb ~ 8 x 10 ~ 6 mag kpc 2 Af^ 1 is the B-band mass absorption 



coefficient (Sadler & Gerhard 1985). The dust extinction in 



the central regions of 3C324 covers an area of at least 2 
by 2kpc, requiring a minimum of ~ lO 6 Af0 of dust to be 
concentrated near the nucleus of the radio galaxy. 

The extended emission of 3C324 is strongly polarised 
due to scattering of radiation from a hidden quasar nucleus, 
and the shape of the polarised flux spectrum indicates that, 
regardless of whether the scattering medium is electrons or 
dust, the scattered light must be dust reddened. In the case 
of pure dust scattering, a mass of at least a few time s 10 6 A/p) 
of dust must be distributed throughout the galaxy ( |Cimatti 
et al. 1996). In addition to scattering the AGN emission, 



this same dust would absorb some of the AGN radiation 
and reprocess this to submillimetre wavelengths. These two 
lower limits allow the mass of warm dust in 3C324 to be 
constrained to within an order of magnitude even if the sub- 
millimetre derived value is considered as an upper limit. 

The star formation rate of 35OM0 yr _1 deduced for 
3C324 from the SCUBA observations can be compared 
with a star-formation rate estimated from the rest-frame 
2800A flux density. The Keck spectrum of Cimatti et al. 
(199m) provides a 2800A continuum flux density of 2.6 x 

in — J.o „ — 1„™ — 2 X — 1 /v.„, — „ o o u,. i ^1^*- 



et al. 1996). The mass of this centrally concentrated dust 



10~ iB ergs _1 cm _ ^A~ i from a 3.8 by 1 arcsecond slit along 
the galaxy. Extinction by a few times lO 6 Af0 of dust spread 
evenly throughout this region would mean that the intrin- 
sic flux would be about 50% higher. However, the strong 
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UV/optical polarisation of this galaxy suggests that between 
30 and 50% of the 2800A flux density of 3C324 is associated 
with a scattered component (Cimatti et al. 1996), and the 
nebular continuum contribution will also be highly signifi- 
cant, so only a fraction < / < 0.3 will be associated with 
young stars. Comparison with the results of Madau et al. 
( [1998D for the 2800A lu minosi ty expected for on-going star 
formation (for a Scalo (1986) IMF with upper and lower 



cut-offs of 125 and 0.1 Mq respectively, and solar metallic- 
ity), indicates that the rate of on-going star formation must 
be below / x 13OM0 yr" 1 . If, on the other hand, the star 
formation were associated with only a short burst induced 
by the r adio s ource, about 5 x 10 6 yr old, the Bruzual & 
Chariot ( 1993 ) stellar synthesis codes show that the mea- 
sured 2800A flux density corresponds to a star formation 
rate of / x 1O5OM0 yr -1 , just consistent with the value de- 
rived from the current submillimetre observations if / lies 
at the upper end of its allowed range. 

With these star formation rate limits, powerful radio 
galaxies at redshift one cannot be undergoing anything more 
extreme than short-lived (< 10 7 yr) starbursts of a few hun- 
dred solar masses per year. This is in stark contrast to the 
situation in radio galaxies at redshifts z > 3 where star- 
formation rates of several thousands of solar masses are de- 
rived. These values are based upon the assumption that the 
dust is heated by young stars rather than by the AGN. The 
current observations support that hypothesis: the AGN of 
the high redshift galaxies have a similar intrinsic power to 
that of 3C324, whilst the submillimetre emission from 3C324 
at least a factor of a few lower. Therefore, only a small frac- 
tion of the dust in the very distant galaxies can be heated 
by the AGN. 

Note added in manuscript: On August 19 1998, a deeper 
(64 sub-scans) Effelsberg 32 GHz observation of 3C324 was 
made during better weather conditions, using the same ob- 
servational set-up as described in the text. An improved 
flux density measurement of 38.5 ± 3.6 mjy was determined, 
strengthening the conclusion of a turn-over in the syn- 
chrotron spectrum. 
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